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Abstract 

Artificial intelligence needs methods for object modeling 
having the set of ley pointe. A novel method of Hurwitz- 
Radon Matrices (MHR) is used in signature modeling. 
Proposed method is based on the family of Hurwitz-Radon 
(HR) matrices which possess columns composed of 
orthogonal vectors. Two-dimensional curve as signature is 
modeling via different functions as probability distribution 
functions: polynomial, sinus, cosine, tangent, cotangent, 
logarithm, exponent, arc sin, arc cos, arc tan and power 
function. It is shown how to build the orthogonal matrix 
OHR operator and how to use it in a process of signature 
modeling. 
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Introduction 

Signature description by key points is a critical part in 
many applications of image processing. Thus artificial 
intelligence have a problem: how to model the curve 
[1,2] via discrete set of two-dimensional points? Also 
the matter of signature representation and description 
is still opened [3,4]. The author wants to approach a 
problem of signature representation by characteristic 
points. Proposed method relies on functional 
modeling of curve points situated between the basic 
set of the nodes. The functions that are used in 
calculations represent whole family of elementary 
functions: polynomials, trigonometric, cyclometric, 
logarithmic, exponential and power function. These 
functions are treated as probability distribution 
functions in the range [0;1]. Nowadays methods apply 
mainly polynomial functions, for example Bernstein 
polynomials in Bezier curves, splines and NURBS [5]. 
Numerical methods for data interpolation are based 
on polynomial or trigonometric functions, for example 
Lagrange, Newton, Aitken and Hermite methods. 
These methods have some weak sides [6] and are not 
sufficient for object modeling in the situations when 



the signature cannot be build by polynomials or 
trigonometric functions. Also trigonometric basis 
functions in Fourier Series Shape Models are not 
appropriate for describing all shapes. Model-based 
vision such as Active Contour Models (called Snakes) 
or Active Shape Models use the training sets to fit the 
data and they are applied only for closed curves. In 
this paper discussed approach is not limited to closed 
curves and it does not use a training set of some 
images, but only a set of two-dimensional nodes of the 
signature. Proposed signature modeling is the 
functional modeling via any elementary functions and 
it helps us to fit the curve and to match the signature 
in object modeling or image analysis. 

The author presents novel method of flexible 
signature modeling. This paper takes up new method 
of two-dimensional curve modeling via using a family 
of Hurwitz-Radon matrices. The method of Hurwitz- 
Radon Matrices (MHR) requires minimal assumptions 
about object. The only information about shape or 
curve is the set of at least five nodes. Proposed method 
of Hurwitz-Radon Matrices (MHR) is applied in curve 
modeling via different coefficients: polynomial, 
sinusoidal, cosinusoidal, tangent, cotangent, 
logarithmic, exponential, arc sin, arc cos, arc tan or 
power. Function for coefficient calculations is chosen 
individually at each signature modeling and it 
represents probability distribution function. MHR 
method uses two-dimensional vectors (x,y) for curve 
modeling. Shape of the object is represented by 
succeeding points (xi,yi) e R 2 as follows in MHR 
method: 

• MHR version with no matrices (N= 1) needs 2 
nodes or more; 

• At least five nodes (xi,yi), (xiyi), (x3,j/3), (x4,j/4) 
and (X5,ys) if MHR method is implemented 
with matrices of dimension N = 2; 

• For better modeling nodes ought to be settled 
at key points of the curve, for example local 
minimum or maximum and at least one point 
between two successive local extrema. 
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Condition 2 is connected with important features of 
MHR method: MHR version with matrices of 
dimension N = 2 (MHR2) requires at least five nodes, 
MHR version with matrices of dimension N = 4 
(MHR4) needs at least nine nodes and MHR version 
with matrices of dimension N = 8 (MHR8) requires at 
least 17 nodes. Condition 3 means for example the 
highest point of the object in a particular orientation, 
convexity changing or curvature extrema. So this 
paper wants to answer the question: how to model the 
signature for discrete set of points? 



FIG. 1 NODES OF THE SIGNATURE BEFORE MODELING 

Coefficients for signature modeling are computed via 
individual features of the object using probability 
distribution functions: polynomials, power functions, 
sinus, cosine, tangent, cotangent, logarithm, exponent 
or arc sin, arc cos, arc tan. 

Probabilistic Modeling 

The method of Hurwitz - Radon Matrices (MHR), 
described in this paper, is computing points between 
two successive nodes of the curve. Signature data are 
interpolated and parameterized for real number 
ae[0;l] in the range of two successive nodes. MHR 
calculations are introduced with square matrices of 
dimension N = 1, 2, 4 or 8. Matrices At, i = 1,2. ..m 
satisfying 

AjAk + AkAj= 0, A? = -I for * k; j, k=l,2...m 

are called a family of Hurwitz - Radon matrices, discussed 
by Adolf Hurwitz and Johann Radon separately in 
1923. A family of Hurwitz - Radon (HR) matrices [7] 
are skew-symmetric (A; T = -At), Ar 1 = - Ai and only for 
dimension N = 1, 2, 4 or 8 the family of HR matrices 
consists of N - 1 matrices. For N = 1 there is no matrices 
but only calculations with real numbers. For N = 2: 
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For N = 8 we have seven HR matrices with elements 0, 
±1. So far HR matrices have found applications in 
Space- Time Block Coding (STBC) [8] and orthogonal 
design [9], in signal processing [10] and Hamiltonian 
Neural Nets [11]. 

How coordinates of curve points are applied in 
Active Object Modeling? If boundary points have the 
set of following nodes {(xi,y;), i = 1, 2, n] then HR 
matrices combined with the identity matrix In are used 
to build the orthogonal Hurwitz - Radon Operator 
(OHR). For point pi=(xi,yi) OHR of dimension N = 1 is 
represented by matrix (real number) Mi: 

x 1 X, 

For points pi=(xi,yi) and f>2=(x2,y2) OHR of dimension 
N=2 is build via matrix M2: 
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For points pi=(xi,yi), p2=(x2,y2), p3=(x3,y3) and pi={x<i,yi) 
OHR Mi of dimension N = 4 is introduced: 
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where 
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For nodes pi=(xi,yi), p2=(x2,y2),... and ps=(xs,ys) OHR Ms 
of dimension N = 8 is constructed [12] similarly as (1) 
and (2): 
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and u = (ho, hi,..., hz) t (4). OHR operators Mn (0)-(3) 
satisfy the condition of interpolation 

MfX = y (6) 
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for x = (xi,X2...,xn) t e R N , x # 0, y = (yi,y2...,yw) T e R N 
and N= 1,2,4 or 8. 

Distribution Functions in MHR Probabilistic 
Mode ling 

Coordinates of points settled between the nodes are 
computed [13] using described MHR method [14]. 
Each real number c e [a;b] is calculated by a convex 
combination c = a ■ a + (1 - a) ■ b for 

a = — £[0;l]. (7) 
b — a 

The average OHR operator M of dimension N= 1, 2, 4 
or 8 is build: 

M = y ■ A + (l-y)- B ■ (8) 

The OHR matrix A is constructed (l)-(3) by every 
second point pi=(xi=fl,yi), p3=(x3,y3),...and piN-\=(xiN-\, 
yw-i): 

A = M.N(pi, pi,..., p2N-l). 

The OHR matrix B is computed (l)-(3) by data 
p2=(x2=b,y2), p4=(x4,y4),... and j)m=(x2N,y2N): 

B = Mn(P2, fi,..., plN). 

Vector of first coordinates C is defined for 

a= a-X2>-i+ (l-a)-X2> , i=l,l,...,N (9) 

and C = [a, a,..., cn] t . The formula to calculate second 
coordinates y(c;) is similar to the interpolation formula 
(5): 

Y(C) = M-C (10) 

where Y(C) = [y(ci), y(c2),..., y(cN)] T . So modeled value 
of y(ci) from (9) depends on two, four, eight or sixteen 
(2N) successive nodes. For example N = 1 results in 
computations without matrices: 

A = M l {p l ) = ^>B = M l {p 2 ) = ^> C = ci= a-xi+ (l-a)-X2, 

r(C) = y( Cl ) = (r A +a-r)— 

y(c,) = a ■ y ■ y t + (l- a)(l- y)y 2 + y(l - a)^-x 2 + a(l - y)^-x l 

(11) 

Key question is dealing with coefficient y in (7). Basic 
MHR version means y = a and then (10): 

y(q) = a 2 -y, +(l-a) 2 y 2 +a(l-a)(^-x 2 +^ 2 -x 1 ). ( 12 ) 

Xj x 2 

Formula (11) differs from linear interpolation: 

y(c) = a-y 1 +(l-a)y 2 . 



www.seipub.org/gpai/ 
MHR is not a linear interpolation. 

Each individual signature requires specific 
distribution of parameter y (7) and y depends on 
parameter a e [0;1]: 

Y = F(a), F(0) = 0, F(l) = l 
and F is strictly monotonic. 

Coefficient y is calculated using different functions 
(polynomials, power functions, sinus, cosine, tangent, 
cotangent, logarithm, exponent, arc sin, arc cos, arc tan) 
and choice of function is connected with initial 
requirements and signature specifications during 
fitting and matching of the object. Coefficients y and a 
(6) are strongly related: 

• y = 0<-+a=0; 

• y = l<-+a =1; 

• ye[0;l]. 

Efifferent values of coefficient y are connected with 
applied functions F(a). These functions (12)-(38) 
represent the probability distribution functions for 
random variable a e [0;1] and real number s > 0: 
power function 

y = a s with s>0. (13) 
For s = 1: basic version of MHR method when y = a. 
sinus 



y = sin(a s ■ nil) , s > 


(14) 


or 




Y = sin s (a ■ n/1) , s > 0. 


(15) 


For s = 1: y = sin(a ■ n/1). 
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cosine 




Y = l-cos(a s ■ n/1) , s > 
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or 




Y = l-cos s (a ■ n/1) , s > 0. 


(18) 


Fors = l: y = l-cos(a ■ n/1). 
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tangent 




Y = tan(a s ■ ti/4) , s > 
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or 




Y = tan s (a ■ ti/4) , s > 0. 


(21) 


For s = 1: y = tan(a ■ ti/4). 


(22) 


logarithm 




y = log2(a s + 1) , s>0 


(23) 


or 




Y =log2 s (a + 1), s >0. 


(24) 


Fors = l: y = log2(a + l). 


(25) 



exponent 

Y = ( a ll y , s > and a > and a * 1. (26) 
a-l 
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y = 2"- 1. 



For s = 1 and a = 2: 
arcsine 
y = 2/n- arcsin(a s ) , s > 
or 

Y = (2/n- arcsin a) s , s > 0. 

For s = 1: y = 2/n- arcsin(a). 

arccosine 

Y = 1-2/71- arccos(a s ) , s > 

or 

Y = 1-(2/tt- arccos a) s , s > 0. 

For s = 1: y = 1-2/n- arccos(a). 

arctangent 

Y = 4/71- arctan(a s ) , s > 

or 

Y = (4/n- arctan a) $ , s> 0. 

Fors = l: y = 4/n- arctan(a). 

cotangent 

Y = ctg(n/2 - a s ■ 7t/4) , s > 
or 

Y = cfg s (ti/2 - a ■ 7i/4), s > 0. 
For s = 1: y = ctg(n/2 - a ■ 7t/4). 



(27) 

(28) 

(29) 
(30) 

(31) 

(32) 
(33) 
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(35) 
(36) 

(37) 

(38) 
(39) 



Functions used in y calculations (12)-(38) are strictly 
monotonic for random variable a e [0;1] as y = F(a) - 
probability distribution function. Choice of function 
and value s depends on object specifications and 
individual requirements. Fixing of unknown 
coordinates for curve points using (6)-(9) is called by 
author the method of Hurwitz - Radon Matrices (MHR) 
[15] . So here are five steps of MHR modeling: 

• Choice of nodes at key points. 

• Fixing the dimension of matrices N = 1, 2, 4 or 8: 
N = 1 is the most universal for calculations (it 
needs only two successive nodes to compute 
unknown point between them) and it has the 
lowest computational costs (10); MHR with N = 2 
uses four successive nodes to compute unknown 
coordinate; MHR version for N = 4 applies eight 
successive nodes to get unknown point and MHR 
with N = 8 needs sixteen successive nodes to 
calculate unknown coordinate (it has the biggest 
computational costs). 

• Choice of distribution y = F(a): basic distribution 
for y = a. 

• Determining values of a: a = 0.1, 0.2... 0.9 or 0.01, 
0.02. ..0.99 or others. 

• The computations (9). 

These five steps can be treated as the algorithm of 
MHR method of modeling and interpolation (6)-(9). 



Considering nowadays used probability distribution 
functions for random variable ae[0;l] - one 
distribution is dealing with the range [0;1]: beta 
distribution. Probability density function /for random 
variable a e [0;1] is: 



f(a) = c-a s -(\-aY ,s>0, r>0 



(40) 



When r = probability density function (39) represents 
f(a) = c-a 5 an d then probability distribution function 
F is like (12), for example f(a) = 3a 2 an d Y = a3 - If s 
and r are positive integer numbers then y is the 
polynomial, for example /(«) = 6a(l-a) an d y = 3a 2 - 
2a 3 . So beta distribution gives us coefficient y in (7) as 
polynomial because of interdependence between 
probability density /and distribution F functions: 

f(a) = F'(a), F(a) = ]f(t)df ( 41 ) 



For example (40): f(a) = a-e a and 
y = F(ct) = (a-\)e a +1- 

What is very important: two signatures may have the 
same set of nodes but different N or y results in 
different signatures (Fig.6-13). Here are some 
applications of MHR method with basic version (y=a): 
MHR2 is MHR version with matrices of dimension 
N=2 and MHR4 means MHR version with matrices of 
dimension N= 4. 



Wykres inierpolacji krzywej MHR-2 



5S 0,0 



FIG. 2 FUNCTION F(X) = X?+X 2 -X+l WITH 396 INTERPOLATED 
POINTS US ING BAS IC MHR2 WITH 5 NODES 



Wykres Intarpolacji krzywej MHR-2 



>• -1 .* 



FIG. 3 FUNCTION F(X) = X 3 +LN(7-X) WITH 396 INTERPOLATED 
POINTS LB ING BASIC MHR2 WITH 5 NODES 
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FIG. 4 FUNCTION F(X) = XH2X-1 WITH 792 INTERPOLATED 
POINTS USING BASIC MHR4 WITH 9 NODES 



Wvkrea interpolacj i krcywej MHR-2 

; . wii » X? — | 



o.a 1.0 1.1 1.2 1.3 M i.s i.e 1.7 i,a i.b 
wsp6irzedfte x 



FIG. 5 FUNCTION F(X) = 3-2 x WITH 396 INTERPOLATED POINTS 
USING BAS IC MHR2 WITH 5 NOD ES 

Figures 2-5 show modeling of continues data 
connected with determined functions. So these 
functions are modeled. Without knowledge about 
function, signature modeling has to implement the 
coefficients y (12)-(40), but MHR is not limited only to 
these coefficients. Each strictly monotonic function F 
between points (0;0) and (1;1) can be used in MHR 
signature modeling. 

Applications of Mhr Probabilistic Modeling 

Signature nodes: (0.1;10), (0.2;5), (0.4;2.5), (1;1) and (2;5) 
from Fig.l are used in some examples of MHR 
signature modeling with different y. Points of the 
signature are calculated for N = 1 and y = a (11) in 
example 1 and with matrices of dimension N=2 in 
examples 2 - 8 for a = 0.1, 0.2,. . .,0.9. 

Signature modeling for N = 1 and y = a. 



3 * 
♦ 

6 * 



FIG. 6 MODELING WITHOUT MATRICES (N = 1) FOR NINE 
RECONSTRUCTED POINTS BETWEEN NODES 



Sinusoidal signature modeling with y = sin(a ■ tl/2). 



FIG. 7 SINUSOIDAL MODELING WITH NINE RECONSTRUCTED 
SENATURE POINTS BETWEEN NODES 

Tangent modeling for y = tan (a ■ tt/4). 



♦ 
♦ 



FIG. 8 TANGENTSIGNATURE MODELING WITH NINE 
INTERPOLATED POINTS BETWEEN NODES 

Tangent modeling with y = tan(a k ■ tt/4) and k = 1.5. 



\ 



FIG. 9 TANGENT MODELING WITH NINE RECOVERED 
SENATURE POINTS BETWEEN NODES 

Tangent signature modeling for y = tan(a k ■ tt/4) and k 
= 1.797. 



* 
♦ 



FIG. 10 TANGENT MODELING WITH NINE RECONSTRUCTED 
POINTS BETWEEN NODES 
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Sinusoidal modeling with y = sin(a k ■ nil) and k = 2.759. Cone lusions 



A. 



U.i 



2,5 



FIG. 11 SINUSOIDAL MODELING WITH NINE INTERPOLATED 
SIGNATURE POINTS BETWEEN NODES 

Power function modeling for y = a k and k = 2.1205. 



V 



FIG. 12 POWER FUNCTION SIGNATURE MODELING WITH 
NINE RECOVERED OBJECT POINTS BETWEEN NODES 

Logarithmic signature modeling with y = log2(a k + 1) 
and/: = 2.533. 



0.5 



1.5 



FIG. 13 LOGARITHMIC MODELING WITH NINE 
RECONSTRUCTED POINTS BETWEEN NODES 

These eight examples demonstrate possibilities of 
signature modeling for key nodes. Reconstructed 
values and interpolated points, calculated by MHR 
method, are applied in the process of curve modeling 
for fitting and matching the signature during its 
analysis. Every individual signature or handwriting, 
each letter or number can be modeled by some 
function for parameter y. This parameter is treated as 
probability distribution function for each individual 
signature 



The method of Hurwitz-Radon Matrices (MHR) 
enables modeling of two-dimensional curves using 
different coefficients y: polynomial, sinusoidal, 
cosinusoidal, tangent, cotangent, logarithmic, 
exponential, arc sin, arc cos, arc tan or power function 
[16]. Function for y calculations is chosen individually 
at each signature modeling and it is treated as 
probability distribution function: y depends on initial 
requirements and curve specifications. MHR method 
leads to signature modeling via discrete set of fixed 
points. So MHR makes possible the combination of 
two important problems: interpolation and modeling. 
Main features of MHR method are: 

• modeling of L points is connected with the 
computational cost of rank O(L) ; 

• MHR is well-conditioned method (orthogonal 
matrices) [17]; 

• coefficient y is crucial in the process of 
signature modeling and it is computed 
individually for each object (letter or figure). 

Future works are going to: features of coefficient y, 
implementation of MHR in object recognition [18], 
shape representation, curve fitting, contour modeling 
and parameterization [19]. 
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